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Spatially uniform electric fields have been used to induce instabilities in liquids and polymers, 
and to orient and deform ordered phases of block-copolymers. Here we discuss the demixing phase 
transition occurring in liquid mixtures when they are subject to spatially nonuniform fields. Above 
the critical value of potential, a phase-separation transition occurs, and two coexisting phases appear 
separated by a sharp interface. Analytical and numerical composition profiles are given, and the 
interface location as a function of charge or voltage is found. The possible influence of demixing on 
the stability of suspensions and on inter-colloid interaction is discussed. 



INTRODUCTION 

Electric fields influence the structure and thermody- 
namic behavior of charged as well as neutral matter. 
Their effect is strong, they can be switched on or off, 
and they are easily scalable to the sub-micron regime 
[l[ . There are two main distinctions with respect to the 
field: spatially uniform vs. nonuniform fields. There 
are also two broad classes of material properties: pure 
dielectric vs. conducting media. All four combinations 
are relevant to phase transitions in liquid and polymer 
mixtures and to liquid- vapor coexistence in pure liquids. 

Perfect dielectrics. The electrostatic energy of di- 
electric materials is given by the expression 
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where e is the dielectric constant and E is the electric 
field. The negative sign before the integral is applicable 
to situations where the electric potential -0 (E = —Vip) 
is given on the bounding surfaces; in cases where the 
charge is prescribed, E is given as a function of the dis- 
placement field D, and the Legendre transform reverses 
the sign 

The phase-transition described below occurs in sys- 
tems described by bistable free-energy functionals 
giving rise to a phase-diagram in the composition- 
temperature plane divided into two regions: homo- 
geneous mixture and a phase-separated state. For 
concreteness, we consider the symmetric mixture free- 
energy density /„ given by 
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This free-energy is given in terms of the dimcnsionless 
composition (0 < < 1). In a binary mixture of two 
liquids 1 and 2, with dielectric constants ei and £2, is 
the relative composition of (say) liquid 2. In A/B poly- 
mer blends, it is the relative volume fraction of polymer 
A, and similarly for an A/B diblock-copolymer melt. 
vq is a molecular volume, ks is the Boltzmann con- 
stant and Tc is the critical temperature. In symmetric 
mixtures, the transition (binodal) temperature Tt at a 



given composition is given by dfm{Tt, 4>) /dip = 0, that is, 
Tt/Tc = 2(20 -l)/log(0/(l-0)) H. In the absence of 
electric field, the mixture is homogeneous if T > Tj, and 
unstable otherwise. The field-induced phase-transition 
discussed below does not depend on the exact form of 
fm\ it occurs also in a Landau series expansion of Eq. 
([2]) around the critical composition 0c, and in other 
forms having a "double- well" shape. 

The dielectric constant e depends on <f> via a consti- 
tutive relation. A variation of from its critical value, 
0c, induces a variation of e from the critical permittiv- 
ity Sc. When the composition deviation ip = (f> — (pc 
is small enough, \ip\ ^ 1, the constitutive relation e(0) 
can be written as a Taylor series expansion to quadratic 
order: 
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The "dielectric contrast" Ae is simply equal to £2 — £1, 
if e" vanishes. 

The electric field depends on the imposed external po- 
tentials or charges and on the local dielectric constant. 
Let us denote by Eq the electric field corresponding to 
the system with uniform composition 0c everywhere. 
Composition changes in induce changes in e, and since 
e and E are coupled via Laplace's equation V(eE) = 0, 
one has variations in electric field. 

We may thus write to quadratic order in ip 
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Note that Eq is constant in space only inside a parallel- 
plate capacitor, or if the sources of the field are very 
far from the system under investigation. Clearly, even 
if Eq is uniform, composition variations lead to field 
nonuniformities . 

One can expand the electrostatic energy density in 
Eq. ([1]) in powers of ip: 
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The unimportant constant corresponds to the electro- 
static energy of the system with uniform composition, 
and it serves as a reference energy. If the field Eq is uni- 
form in space, the two terms in linear order of Lp simply 
add a constant to the chemical potential, and therefore 
are inconsequential for the thermodynamic state of the 
system 

Landau and Lifshitz showed that the existence of a 
e"ip'^ term in Eq. (O is responsible to a shift of the 
critical temperature Tc [H, [J. They found that is 
increased by ATc given by 
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Tc and the whole binodal curve close to the critical point 
are increased if e" > (field-induced demixing) or de- 
creased if e" < (field- induced mixing). A similar ex- 
pression exists for a pure liquid in coexistence with its 
vapor. 

The experiments, starting with P. Debye and Kle- 
both [il, are in contradiction with this prediction. De- 
bye and Kleboth investigated the critical temperature 
of a Isooctane-Nitrobenzene mixture (relative permit- 
tivities 2.0 and 34.2, respectively) They observed re- 
duction of Tc by 15 mK in a field of 4.5 V/fim. Their 
measurements were later verified by Orzechowski 6]. 
Beaglehole worked with a Cyclohexane- Aniline mixture 
(relative permittivities 2 and 7.8, respectively), and he 
measured reduction of Tc by as much as 80 mK in a 0.3 
V//im dc field Early worked on the same mixture 
but in 1 V//im ac field, and found no change in Tc- He 
attributed the results of Beaglehole to spurious heating 
. Wirtz and Fuller performed similar experiments on 
n-hexane-Nitroethane mixture (relative permittivities 2 
and 19.7, respectively), and found a reduction of Tc by 
20 mK in a 5 V//im field In all cases, e" was pos- 
itive but still mixing was observed. In addition, the 
observed change in Tc is quite small, typically in the 
10-20 mK range. The only exception is the work of 
Gordon and Reich fu^ . They worked on polymer mix- 
tures of poly(vinyl methyl ether) (PVME)-polystyrene 
(PS) system (relative permittivities 2.15 and 2.6, re- 
spectively), and observed changes significantly larger 
than 1 K. Their strong effect can be attributed to the 
large molecular weight of the polymer (14, 000 - 30, 000 
gr/mol) and to their reduced entropy compared to that 
of simple liquids. 

One is inclined to explain the experimental findings 
by the second and third terms on the second line of Eq. 
([5|) (proportional to (^^). The third term is twice as 
large as the second one and opposite in sign, and the 
two sum to give a free energy contribution proportional 
to the dielectric contrast squared -|-(Ae)^. This is a 
free energy penalty for dielectric interfaces perpendicu- 
lar to the external field. Indeed, these additional terms 
are responsible to the normal field instability in liquids 
and to orientation of ordered phases ( e.g. block- 
copolymcrs) in external fields [HI, 0, E, [3, El, [H, 17]. 
In liquid mixtures they favor mixing (lowering of Tc). 



MIXTURES OF NONPOLAR LIQUIDS IN 
FIELDS GRADIENTS 



Field gradients are general, and occur in all electrodes 
unless special care is taken to eliminate them (super- 
flat and parallel conducting surfaces). When mixtures 
of pure dielectric liquids are subjected to a spatially 
nonuniform field, the situation is very different. The 
direct coupling between field variations and composi- 
tion fluctuations then leads to a dielectrophoretic force, 
depending on Ae in Eq. ([3]), which tends to "suck" the 
component with large e to regions with high electric 





FIG. 1: Three model systems where field gradients lead to 
demixing. (a) A wedge comprised of two fiat electrodes with 
an opening angle /3 and potential difference V. Ri and R2 
are the minimal and maximal values of the distance r from 
the imaginary meeting point, (b) A charged wire with radius 
Ri , or two concentric cylinders with radii Ri and R2 . (c) A 
single charged colloid of radius _Ri and surface charge a. 



Statics 

Three "canonical" geometries with electric field gradi- 
ents are presented in Fig. [T] The first is the "wedge" 
capacitor, made up from two flat and nonpar allel sur- 
faces with potential difference V, and opening angle (3. 
The electric field is then azimuthal, E(r) = V/{(3r)d, 
where r is the distance from the imaginary meeting 
point of the surface, r is bounded by the smallest 
and largest radii Ri and R2, respectively. The sec- 
ond model system consists of a charged wire of radius 
or two concentric metallic cylinders with radii Ri 
and R2 > Ri- In this case the azimuthally-symmetric 
field is E(r) ~ a'Ri/{re{r))r, where a is the charge per 
unit area on the inner cylinder. Lastly, for a charged 
spherical colloid of radius i?i and surface charge tr, 
one readily finds the spherically-symmetric field to be 
E(r) = aRi/{r^e{r))r. 

In all three cases a general scenario occurs: when 
T is above Tc, the composition profile 0(r) is smooth, 
and its gradients increase as the charge on the objects 
increases. However, below Tc the behavior is different - 
0(r) is smooth as long as the charge (voltage) is small, 
and becomes discontinuous when the charge (voltage) 
attains a critical value. At this charge, a sharp interface 
appears between the coexisting domains [T^H^. As the 
charge further increases, the interface location and the 
compositions of the coexisting domains change. 

To see this, consider the wedge capacitor, for which 
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FIG. 2: (a) Graphical solution to Eq. ((Sj. Solid curve is 
/^(<p). Its roots are the transition (binodal) compositions. 
The intersection between f'm{^) and the horizontal dashed 
line gives the solution ip{r) to Eq. (|8]). For voltages V below 
the critical value V* , the dashed line is bounded by lines a 
and b, corresponding to the the maximal and minimal values 
of the right-hand side of Eq. (|8]), giving rise to a continuous 
profile ip(r). At V > V* , line b is displaced to b', and the 
intersection is at < for large r's and at <^ > at small 
r's. (b) Qualitative composition profiles 4>{r). Horizontal 
dashed line is the average composition in the absence of 
field. cjf>(r) varies smoothly when V < V* , and has a sharp 
jump at r = i? when V > V* . 
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the electrostatic energy density is 

Note that we have used a linear constitutive relation. 
In uniform electric fields, such a linear relation would 
mean that the electrostatic energy is simply a constant 
independent of the composition profile. In addition, in 
the wedge geometry the electrostatic energy does not 
have a term proportional to (Ae)^ because the electric 
field is parallel to the dielectric interfaces (both are in 
the 9 direction). 

The equation that governs the composition profile 
<^(r), derived from the Euler-Lagrange equation S{fm + 
fcs)/S(p = 0, is the following: 

2 

+ ^i . (8) 
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fM = -Ae 



Here fi is the chemical potential of the large reservoir at 
infinity. Note that Eq. ([5]) gives an analytical expres- 
sion for r as a function of ip. 



The right hand side of the equation is independent 
of If, and is indicated by the horizontal lines in Fig. [2] 
(a). Suppose the mixture composition in the absence 
of field corresponds to a point above and to the left of 
the binodal curve (homogeneous mixture). A graphical 
solution of the governing equation is obtained by the 
intersection of the horizontal line, whose location de- 
pends on the field, and therefore on r, with the curve 
f^{ip). If T is above Tc, fm is convex, and therefore 
the intersection of the two curves changes smoothly as 
r decreases {E increases). The resulting composition 
profile is shown in Fig. [2] (b) . 

However, the situation is different below Tc- here 
fmi'P) behaves like —ip + (p^. When the applied voltage 
is small enough such that the maximum value of the 
right-hand side of Eq. ([8]) occurs at line b of Fig. [2] 
(a), as one goes from large to small values of r (increas- 
ing E), the composition increases, but ip{r) is always 
continuous. There is a critical value of the voltage, V*, 
where this is not true: above the critical potential, the 
maximum value of the horizontal line can be at b' in the 
figure. Therefore, (p increases with decreasing r until, 
at a certain location r = R, there are three solutions. 
The middle one is an unstable while the other two are 
stable. At this point, the composition "jumps" between 
the two stable values and a discontinuity appears. At 
such voltages, the profiles are discontinuous and the co- 
existence between two distinct phases occurs. 

Assuming that the jump in (/) occurs at the binodal 
values, one obtains the stability criterion [l9| 
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Here E — V/{(3Ri) is the largest value of the field. A 
mixture of initial homogeneous composition (^o is unsta- 
ble and demixes into two coexisting domains under the 
given field if the temperature is below Tt + AT, where 
Tt{4>o) is the zero- field transition (binodal) temperature 
at composition 00 . In contrast to uniform fields, where 
field variations result from composition variations, here 
field gradient are due to the non-flat geometry of elec- 
trodes. Hence, AT above is typically 2-100 times larger 
than AT in uniform fields [Eq. ([6])]. Note that simi- 
lar demixing is also expected to occur in a rapidly ro- 
tating centrifuge. In that case (wr)^ is the analogue of 
the spatially-dependent field E"^ , where r is the distance 
from the rotation axis and w the angular frequency. The 
density difference Ap = p2 — pi replaces the dielectric 
contrast Ae [2ll |. 

Eq. ([9|) may be inverted to give the critical volt- 
age for demixing V* as a function of 0o and tem- 
perature. One finds that oc (T - Ttf"^. In 
the experiments of the Leibler group, conducted us- 
ing sharp "razor-blade" electrodes, the measured ex- 
ponent was 0.7 ± 0.15, larger than the value 1/2 cited 
here. One may write the dimensionless potential as 

1/2 

= V[v„e„/{Ap^kBT,Rl)y'\ where £0 is the vac- 
uum permittivity. The critical value of Uw for a closed 
wedge, U^, is obtained by an approximation similar to 
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that of Eq. 
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i]), namely [22] 
ksT, A\Ae\/eo dc^^ 



(10) 



where 0t is the transition composition, x = R2/R1, 
and the dimensionless function g is given by g{x) = 
2(^2 - l)/(a;2 - 1 - 21na;). 

Dynamics 

The phase ordering dynamics of mixtures in electric 
field is quite different from the no-field case, since the 
electric field introduces a preferred direction and thus 
breaks the initial system symmetry. The phase transi- 
tion studied here is even more difficult, because spatially 
nonuniform fields also break the translational symme- 
try. 

In this phase transition, droplets nucleate every- 
where, not only in regions of high electric field. As they 
grow, they move under the external force. The viscos- 
ity plays an important role, in addition to the field's 
amplitude, location in the phase diagram and distance 
from the binodal, and dielectric constant mismatch 
Ae. Clearly, the spatial dependence of the electric field 
means the initial destabilization and phase ordering dy- 
namics are quite different from the well-studied normal- 
field instabihty in thin liquid films '3, 'is', H, [H, [l!] and 
the regular coarsening dynamics [27. , ,28i] . 

For salt-free mixtures, the starting point for the dy- 
namics is the following set of equations 2^, 2^, 30t : 
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u is the velocity field corresponding to hydrodynamic 
flow and 77 is the liquid viscosity. Equation (|lip is a con- 
tinuity equation for 0, where —L'S/{df/d(j)) is the diffu- 
sive current due to inhomogeneities of the chemical po- 
tential, and L is the transport coefhcient (assumed con- 
stant). Equation is Laplace's equation, Eq. (fT5|) 
implies incompressible flow, and Eq. (jl4p is Navier- 
Stokes equation with a force term —(pS/df/dcj) (27l[28t. 
It should be noted that similar models have been pro- 
posed in the literature; the main differences here are 
the bistability of /„i and the nonuniform fields derivable 
from the potential ip. The presence of salt is naturally 
incorporated into the model by adding two continuity 
equations for the two ionic species, and by using Pois- 
son's equation instead of Laplace's. As a starting point, 
we assume there is no net flow due to pressure gradients 
or moving solid surfaces - flow will be purely a result of 
the forces exerted by the electric field. In addition, the 
liquid viscosity 77 is taken as a simple constant scalar, 
independent of mixture composition. 

Since the equations arc coupled and nonlinear, it is 
useful to first study the demixing in one of the ge- 
ometries mentioned above (sphere, cylinder, or wedge). 



Consider, for example, the simplifications of the phase- 
ordering equations occurring in the system of concen- 
tric cylinders. In this annular capacitor, the no-flow 
conditions on the inner and outer cylinders lead to 
a vanishing flow velocity: u = everywhere. One 
is therefore left with only a single equation to solve, 
84)/ dt — LS/'^Sf / 54). This equation can be viewed as 
a continuity equation d4>/dt = —V • J, with a current 
density J = —LVSf/S4- In a closed system, J vanishes 
at Ri and R2, and the integral 27rr0(r, t)dr is kept 
constant throughout the system dynamics. Gauss's law 
readily gives the electric field in the concentric capac- 
itor when the charge is given. In the explicit scheme 
we used, 4'(t) is given by a successive summation of 
— (V • 3)dt, calculated for each time interval dt. The 
initial condition for the calculation was a homogeneous 
distribution 4>o- 

When the charge on the capacitor is larger than a 
threshold charge, we observe fast creation of a disconti- 
nuity near the inner cylinder which then starts to move 
outsides. The location of this front, separating the inner 
and outer regions, i?(i), is shown against time in Fig. [3] 
(a). In closed systems, R cannot grow indefinitely, since 
mass conservation dictates an upper bound i?max given 

by 
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In the numerical calculation, we find a match to an 
exponential relaxation with a single time constant t: 



R{t) — Rt=o + (Roo — -Rt=o)(l 
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Roc corresponds to the steady-state solution; it tends to 
^max when the voltage or charge tend to infinity. The 
time constant r depends on the external potential (or 
charge) and on the temperature and composition. Part 
(b) plots r as a function of temperature for different 
average compositions. The calculations indicate faster 
dynamics (smaller r) when the average composition is 
farther from the critical value (large \4o ~ 4'c\) or when 
the cylinder's charge is large, see Fig. [3] (b). 



MIXTURES CONTAINING SALT 

Mixtures of polar liquids (e.g. aqueous solutions) 
contain some amount of charge carriers. In such mix- 
tures, the physics is rich and quite different from the 
simple dielectric case. The most important feature is 
due to screening, occurring when dissociated ions ac- 
cumulate at the charged surfaces. This means that 
the electric field is substantial only close to the sur- 
faces, within the screening distance A. Field gradi- 
ents thus originate from both geometry and screen- 
ing, and the phase transition depends on at least two 
lengths. The ionic screening therefore adds to the di- 
electrophoretic force which separates the liquids compo- 
nents from each other. Since screening is omni-present, 
phase-separation may occur even near parallel and fiat 
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FIG. 3: (a) Plot of the R{t), the dynamics of front location 
between coexisting phases, for a mixture confined by two 
concentric cylinders when the charge on the inner cylinder 
is above the demixing threshold. <j)o = 0.3, T/Tc = 0.95, and 
the dimensionless charge of the inner cylinder is Uc = 0.445, 
where Uc = a [vq/ (ieokBTc)]^^^ , a is the charge per unit 
area of the cylinder, and Ri is its radius. The exponential 
time constant r in Eq. (|16|l is plotted in (b) for different 
values of 00 and Uc- Blue circles: 0o — 0.2, Uc = 0.252. 
Green stars: = 0.3, Uc = 0.199. Red squares: 0o = 
0.3, Uc ~ 0.252. R is scaled by and time is scaled by 
LkBTl{v^,R\). 



charged surfaces, i.e. in one dimension. But ions have 
another effect besides increasing the dielcctrophorctic 
force. Ions have in general different solubilities in the 
different liquids. As an ion drifts toward the electrode, 
it might "drag" with it the preferred liquid component 
[31I, [s^. Thus, the solubility introduces a force of elec- 
trophoretic origin, proportional to the ions' charge. 

We use the following free energy density to describe 
the system: 
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+ fcsT [n^ In (uon^) + n In (fon )] — [i<i> 

— (Au^n^ + Au^rt^) — A^n^ — \^ nr + const. 

The free energy depends on three fields: the electric po- 
tential ip{r), and the two number densities of positive 
and negative ions: n(r)*. The new terms added here 
are the interaction of ions with the potential (n^?/') s-nd 
the ideal-gas entropy of ions (logarithmic terms). In 



addition, the parameters Au+ and Au^ measure the 
affinity of the positive and negative ions toward the 
liquid- 1 environment, respectively [33| . Am"*", for exam- 
ple, measures how much a positive ions prefers liquid-2 
environment over that of liquid 1. and /x are the La- 
grange multipliers (chemical potentials) of the positive 
and negative ions and liquid composition, respectively, 
and e is the electron charge. 

The free energy is extremized with respect to the 
fields (j), and n^: 
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in keeping with a fixed mixture and ion concentrations: 

(21) 
(22) 



j n±(r)dV = no 



Here V is the total volume and no the average ion 
concentration. The Poisson-Boltzmann equation is ob- 
tained from substitution of Eq. ((20| in the Poisson 
equation, Eq. (fT9|) . 

Due to these forces, the phase transition is expected 
to be greatly enhanced compared to the no-ions case: 
it should occur at elevated temperature above the bin- 
odal, and lead to a very thin demixing layer around the 
charged object [33| . Consider a mixture in the semi- 
infinite space a; > confined by one wall at 2; = 
charged at potential V . An approximate formula for 
the temperature below which a phase transition occurs, 
Tt -\- AT, can be obtained by performing a first loop in 
a perturbative solution of the equations, namely using 
a uniform dielectric constant in Eqs. (fT9| and ((20|) and 
substitution in Eq. ((T8|) . The expression for AT is then 
found to be 13311: 
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Here Am = |Au^|. In most cases, As/ec ^ Au/ksT ^ 
1, and therefore the dielectrophoretic and solubility 
forces have the same magnitudes. The numerator nowo 
is quite small: if we take vo = 8 x 10~^^ m"^ and av- 
erage ion density no = 6 x 10"'^^ m^'^ (lO^^M) we get 
no^o ~ 5 X 10^^. However, AT is usually large. Even if 
we ignore the denominator \(j)c — (j>c\~^ , the exponential 
factor can be huge: if the surface potential is only 1 Volt 
and Tc is the room temperature, we get eV/kgTc ~ 40. 
This shows us that demixing should be observed even 
if the surface potential V or the charge density n^ are 
much smaller. 
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We would like to note that for homogeneous dielec- 
tric hquids, the equation V(eVV') = means that an 
increase of the potential on the bounding surfaces sim- 
ply increases the potential ip proportionally, but that 
for ion-containing mixtures this is not true: due to the 
nonlinearity of the problem, increase of the external po- 
tential leads to a change in the whole distribution ipi^)- 
The composition difference between coexisting phases 
increases with and the front separating the domains 
may move to larger or smaller radii. 



Acknowledgment 

We thank L. Leibler and F. Tournilhac for help in de- 
veloping the ideas presented in this work. This research 
was supported by the Israel Science foundation (ISF) 
grant no. 284/05, and by the German Israeli Founda- 
tion (GIF) grant no. 2144-1636.10/2006. 



CONCLUSIONS 

The steady-state and dynamics of phase transitions 
due to inhomogeneous electric field are discussed. In 
nonpolar mixtures, the composition profile of a mixture 
is given for three geometries with azimuthal or spherical 
symmetries. Above T^, the profile is smooth, while be- 
low Tc it becomes discontinuous if the surface charge or 
voltage exceed their critical values. The location of the 
front separating the two coexisting phases in equilib- 
rium moves to larger values of r as the charge or poten- 
tial increase. In the restricted case shown here, the main 
feature of the dynamical process towards equilibrium is 
the exponential relaxation of the front location. The 
exponential time constant decreases when the poten- 
tial is diminished or when the distance from the critical 
composition is reduced. 

When salt is present, the phase transition is enhanced 
because ionic screening leads to a dielcctrophoretic 
force. In addition, the ions' solubility leads to a strong 
force of electrophoretic origin. Thus, the transition is 
g strengthened and should occur at virtually all tem- 
peratures above the binodal even at modest salt con- 
tent. Our results nicely complement the recent studies 
by Onuki and co-workers on the solvation effects of ions 
in near critical mixtures 18, 2^, 31, 3^ . 

A similar phase transition was observed for a mono- 
layer of surfactant mixture subject to an electric field 
emanating from a charged wire passing perpendicular 
to the monolayer (20| . The more polar surfactant was 
attracted to the wire when the field was applied, while 
the less polar surfactant was repelled. The effect ob- 
served was linear in electric field because (i) the dipoles 
were fixed and not induced, and (ii) they were confined 
to a plane and could not twist up-sidc-down when the 
field's polarity was reversed. 

We point out that when charged colloids are dispersed 
in aqueous solutions, a thin wetting layer could be 
formed due to field-induced demixing, depending on the 
average salt content, temperature, and colloid charge. 
According to Eq. (p3|) . this demixing is quite favorable, 
and one needs not be very close to the binodal curve. 
This should have implications on colloidal aggregation 
[s^ l and on the interaction between charged surfaces in 
solution 'sB'l , because the electrostatically-induced cap- 
illary interaction between the surfaces is expected to be 
attractive. 
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